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I FOREWORD 

As part of its general research program, the Radio Direction Find~ 
ing Redeareh Laboratory at the University of Illinois has undertaken to 
study new systems for radio direction finding. One of the new systems 
proposed. was a "Doppler—Effect" system, The initial theoretical work by 
this group on this system was carried out by Mr. J. Lionel, Boulet and . 
was written up in the Masters! Degree ee opera to in this report. 
Unfortunately, due to a prolonged illness, it became necessary for 
ur. Boulet to withdraw from school and further work on the problem, The 
work was continued by Mr. John Anderson and Mr. Thomas O'Meara, and this 
present report has been written by them. 

Soon after the problem was undertaken, it was found that some experi— 
mental work on this idea had already been carried out at Camp Cole Signal 
Laboratory, the work being summarized in a series of patent applications 


by Paul Hansell. Considerably later, there came to the attention of this 


(3) 
group the British report by C. W. Earp and R. M. Godfrey on "Radio 


Direction Finding by the Cyclical Differential Measurement of Phase." | 
This excellent paper discussed in some detail certain aspects of the 
theory and described a working system which had been built. 

The theoretical approach in this report is aimed chiefly at determin~ 
ing the ability of the system to discriminate against interfering waves 3 
such as those produced by site reflections. Several definite conclusions. | 
in this respect have been drawn and are stated herein. A 'model' Doppler 
Effect system has been built and is in operation. At present , it is 
peing tried out and compared with a simple Adcock system on the Labora~ 
tory's R-D-F System Analyzer. The results of these tests are to be given 


in a later report. 
EK. Ce Jordan 
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Il ABSTRACT 

This report develops and discusses the basic theory of the Doppler 
Type of Direction Finder, The most practical antenna system for this 
purpose appears to be a circular array of fixed antennas with gradual 
or abrupt switching from antenna to antenna around the circle. 

A mathematical discussion of site errors is included for the 
purpose of showing the effect of antenna aperture upon reduction of 
such errors, It is shown that the Doppler Direction Finder offers no 
advantage over the Adcock System when the antenna aperture is a small 
fraction of a wavelength, but that an increase in aperture of the 
Doppler soon reduces site errors to the order of instrumentation errors. 
This comparison with Adcock Systems is made on a figure of merit basis. 

Some discussion of:antenna commutation difficulties and number of 
antennas necessary to place the Doppler in the class of wide aperture 


systems is included. 
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To correlate this with the previous discussion it must be | 
recalled that the rate of change of phase with time is proportional to 


the instantaneous frequency. That is, 


1 do ‘ 
fist ct On dt (8) 





where § for this case = uj t + cos ( 6 - u,t) 


Then 
Taw. 
20 fie = We + ak sin (@- ut ) (9) 
nag, : 
Pee = t3 + ae sin (8 = ult) (10) 


The instantaneous frequency may be expressed as a function of 
the linear velocity, v, =, r of the antenna and as a function of the 


free space propagation velocity,c = Af, , of the incoming wave. 


fa ctl — sin (6 ~ u,t)] (2) 


Notice the similarity between (11) and (5), the only difference 
being that the relative velocity is varying as a sinusoidal function. 


The expression (7), if phase detected, will yield a term 


d 
— cos (8 - wt ) (12) 


which is sinusoidal in nature at the frequency of rotation of the antenna, 
and has a phase dependent upon 9, the angle of arrival of the incoming 
wave. It is possible to obtajn a sinusoidal voltage from the antenna 


motion directly 


nd 
ae cos (wt ) (13) 
which when compared with (12) will give a means of determining ©, Thus 


the basic requirements for a direction finder are fulfilled, 
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(b) Non-Cooperative Systems 
In order for the direction finder to be classified as a Doppler 
type it must possess the equivalent of a moving antenna element of some 
form. However, this does not stipulate the path the antenna must follow. 
® In the thesis, previously mentioned, Mr. Boulet analyzed the case of an 
antenna moving at a uniform speed on a linear path and retracing its 
course periodically. This produces a frequency change depending upon 
speed of motion and upon the angle between the path of the antennae and 
the incoming wave. If the speed is constant the frequency change is a 
measure of the angle. One difficulty arises when the signal arrives 
from the same angle on the other side of the axis. To remove this 
ambiguity a second antenna is made to traverse a path perpendicular to 
the first. 
What appears to be the most practical case to consider is the 
: revolution of the antenna about a vertical axis as was discussed in the 
introduction. Greater linear circumferential speeds are obtained in 
| this case and the problem of getting the antennae back to their starting 
: points in infinitesimal time does not enter, as it does in the linear 
path. Certain transients present in the linear case are not present in 
the circular case, making for smoother operation. 
In the final analysis the rotating antenna will be replaced by 
a circular array of antennas, with electronic switching from antenna to 
antenna. 
! (c) Cooperative Systems 
These two types of Doppler Direction Finders are classified as 


\ 

{ 

: non-cooperative, since no knowledge is known of the nature of the trans- 
| 

mitting source. It should be mentioned that this type of direction find- 
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er can be converted to the so-called beacon system by reversing the role 
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of the rotating antenna from receiving to transmitting. It becories nec 


essary to transmit a synchronizing pulse, which will be received along 
with the phase modulated wave. Since the phase of the phase envelope 
will vary according to the azimuthal angle from the station, a comparison 


with the syrichronizing pulse will give the radial line from the station, 
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LV CIRCULAR ANTENNA ARRAY FOR DOPPLER 
(a) Advantages 
Consider the aforementioned rotating antenna and pick a set of 
representative values for the variables involved. Assuming a rotating 
arm, one meter in length, an angular velocity of 60n radians/sec, 
(1800 R.P.M.) the maximum instantaneous frequency difference will be 
given by (9) 


Af = ,628 x 107° £. cycles/sec. (13) 


max Cc 


The immediate important observation is the small value of Af | 
unless the frequency involved is very large, or the practical Length of 
the arm is exceeded. 

The possibility of producing an effect similar to that described 
for the single rotating antenna, by using fixed antennae connected 


successively to a receiver was investigated. This system should be strict- 


ly referred to as indirect Doppler since the motion is simulated by rapid 
switching. One of the advantages of this indirect approach is the speed 
at which electronic switching can be accomplished, thus eliminating the 
problem of small Af. In addition, the use of fixed antennas makes it 
possible to build arrays of large aperture—a very necessary requirement 
for accurate bearings in the presence of site reflections. Both British . 
and American systems to date utilize a circular array of antennae with 
electronic switching. To appreciate the shape of phase variation involv- 
ed, plate (3) shows an array with a signal arriving at an angle of 45° 
azimuthal, and also the phase of the r.f. voltage picked up as commutation 
progresses around the circle. If the number of antennae used is sufficient, 
the step function of phase, when detected, will yield a fundamental sinu- 
soidal variation similar to the rotating antenna. This variation will 


change phase with the angle of incoming wave in a linear manner, again giv~- 
ing the basic requirement for a Doppler System. 
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(b) Law of Coupling 


The abrupt switching from one antenna to another presents problems 
in detection inasmuch as discrete jumps in phase cause accompanying in~ 
finite frequency change, Eqs (8). To alleviate this a law of coupling to 
the antennae was studied. Consider two elements fn a circular array of 
fixed antennae. Call one element the "nth" element and the next, the 
"(in + 1$h) " element. It is desired to find a law of coupling as switching 
is completed between the two antennae such that the voltage output will 
have a sinusoidal variation of phase, let the voltage received by the 


'tnth't antenna be 


On ( 


eae and (14) 





e, = cos { wtt 
by the "(n+ 1)th" antenna be 


2nn 


nr 








en+] 7 cos [ wt + cos (6 - )] (15) 


where the notation is the same as plate (3). 
The addition of these two voltages with the appropriate law of 
coupling would be required to take the form: 














A ape [ wt + cos ( 6 — 2nat)] (16) 
That is: 
Qnr | 2n(n-1 i 
G(t) cos{wt + — cos[ 9 - oath} } (17) 
Tr 2itn 2nr 
+ H(t) cos[wt + x cos (6 - or) = A cos{wt + 7 cos (6 - 2nat)] 


Consider plate (3). The difference in phase between e, and e, 4) is 


p= af coslO - ty - cos{9 - a); 2 (18) 
and the phase angle (8) of the resultant 
R is given by 
sin = ut) =) sine (19) 
ae 7 
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H(t) e,t Sin 9 





: 2 20 
nan? G(t) e, + H(t) 6,4, cosp (20) 
’ If the arbitrary phase of the 2, voltage is called y the phase 
of R at any instant t, where eee Re is giv en by 
‘ ka ka : 
i H(t) sin 9 
o = + —_—_ 
te) Sy tare cant aye + HE) amo (21) 
Cntt 
ti = constant =c 
Cntt 
This variation must be equal or close to the variation 
a cos(6 ~2nat).The integral to be minimized is then: 
J (a cos(® - 2Qnat) - &(t)] dt 22 
oh {a cos ( at) (t)J (22) 


ka 


The necessary condition for a minimum must satisfy Euler equa~ 


_ tions, that is if in Glx,v,y') ax is to be minimized 
. d OG _ d , 96, _ 2G (23) 
dx @y' y OE see eee Sy" Ox 


Must be satisfied, 


| In our case t=», G(t) 
eS 2 ny es (24) 
| H(t) 
i - sing (25) 
i Ge a) es sin. 25 
| a cos ( Tat) - y ~ are tan [ Seas } 

then ee ye eee (26) 


dt OF '(t) OF (t) 


Since F! (t) is not written explicitely in G 
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Because c 7 0,sing = 0 which leaves F (t) undetermined. Hence, there 


appears to be no general law of coupling which will give the desired 
result. However, Boulet has shown (4) that linear variations of G (t+) 
and H (t) will give a phase curve which closely approximates the ideal 


curves 
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V_ pISCUSSION OF SITH ERRORS 


(a) Constant Phase Front Surface 


The errors encountered in Direction Finders can be classified 


in general as follows: 


(1) Propagation Errors 





The greatest source of this type of error is the 
nonuniform nature of the ionosphere which gives 
rise to a lateral deviation of the main or de~ 
sired ray» For the most part these errors are 
random and may vary rapidly (within the minute) 
or slowly (within the hour). The approach to 


the elimination of errors due to these variations 





seems to be the application of mathematics of 
probability. Spot readings are taken over a 
‘ period of time and the reading of greatest 
probability calculated. The ionosphere can also 
be responsible for the arrival of several rays 


of slightly different azimuth. Minimizing of 


{ 
{ 
| 
| 
| 
| 
| 
| 


these errors also follows statistical lines, 
(2) Instrumentation Errors 

These are errors that occur in the equipment 

and as a result of the type of equipment used. 

Polarization errors, that is errors due to random 

polarization of the incoming signal, are in 

general classified under this heading since their 
: effect can be virtually eliminated by proper 

antenna arrays. Errors due to observation by 


the operator can also fall under this heading. 


med Ores 








(3) Site Errors 
These are errors due to reradiation and reflection 
from objects in the localeof the direction finder. 
Scattering and reradiation by ionic clouds are 
classified here since their effect is the same as 
that for near-by objects. This section is mainly 
interested in the subject of site errors since 
their effect can be materially reduced by the 
choice of the proper type of direction finder. 

When examining the field due to a single incoming wave, a con 
siderable distance from its source, it will be found that the surface 
defining a constant phase is a plane perpendicular to the direction of 
travel, Any direction finder that operates on the phase front principle, 
(such as the Doppler) will indicate an errorless bearing in the direction 
from which the wave is arriving. However, if there is present another 
Signal of appreciable amplitude and of arbitrary phase and direction of 
arrival, the phase front indicated above will no longer be plane. The 
effect of the interfering signal is to corrugate the surface in » manner 
depending upon the relative characteristics of the two signals. When 
the resultant voltage due to the space addition of the two signals is 
examined, it is found to describe a standing wave pattern with points of 
minimum and maximum amplitude spaced at uniform intervals. The phase 
of the resultant voltage also follows definite space patterns and, to 
get. a line of constant phase, it is necessary to set the expression for 
phase at any point equal to a constant and plot the resultant Surves 
This procedure has been outlined in a peesious technical eer cee 


by this laboratory and only the resultant expression is given here: 
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where n = {Q, ly 25 3 ecececes ~l, Ly ~3 waveee) 


ny 
i 


ratio of amplitude of signal #2 to signal #1 


§ 


and 28 <= azimuthal separation of the two rays. 

A plot of constant phase surfaces is shown on plate (5), this 
plate being a reproduction of plate #1 in the above mentioned Summary 
Technical Report. 

(ob) Similarity of all Small Aperture Systems 

If a direction finder with smail aperture, that is small 
dimensions of antenna array, were operating at different points along 
the constant phase front surface, it is clear that the direction of 
-arrival indicated would be a perpendicular to the tangent of the phase 
front curve at the specified point. Since this phase front curve is 
not a straight line, different points would indicate different direc~ 
tions of arrival, It will be seen in a later section that maximum 
error in bearing will occur at points of maximum and minimum amplitude 
of standing wave pattern and that at one point between a maximum and 
minimum there will be indicated the true bearing. From phase front 
considerations, it can be concluded that any direction finder operat~ 
ing on phase front principle and with aperture small in wave length 
will be subject to the same site errors. It might seem that a direc- 
tion finder using an antenna system with a very narrow beam might be 
an exception to this statement if the antenna could be made to have 
small aperture, However, it has been shown* that such antennae 
systems are not practical, because of extremely low sensitivity. Thus 
to get the desired sensitivity, it must be made large in wavelengths 
and falls in the group of wide aperture systems. 

* "Small High Gain Arrays for Direction Finding", N. Yaru, Technical 
Report No.6, Univ. of Ill. Direction Finding Research Laboratory, 
September 1, 1948. 
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Now if the aperture of the system is made larger such that the 


antenna operates over a greater distance on the phase front curve, the 


| 
| 
| 
maximum error variations would be less, This is somewhat analogous to the | 
length of wheelbase on an automobile, Ona road of given roughness, the 
longer the wheelbase, the smoother the ride. The Adcock System does not : 
lend itself to wide aperture operation due to ambiguities which arise. | 
However, the Doppler System is theoretically not restricted on this | 
account and can be increased in aperture to the limit of physical struc~— 
ture and electronic equipment techniques. This tends to put it in the 
class of wide aperture systems of which there are several other examples. 
The site errors of the Doppler System under different apertures will be 
fiven in a later section. 

When increasing the aperture of the Doppler, from a practical 


| viewpoint, the two problems encountered when fixed antennas are used ares 


where to sample the field, and how to detect the resultant phase variation. 





Consider the phase front diagram (Plate 5) and notice that with a Doppler 





z System of only four antennas, properly spaced, it is possible to get no 
phase variation whatsoever as the receiver is switched from antenna to 
antenna. That is, locate the antennas at points of equal phase. This 

| is an extreme example, but it is evident that to get a phase variation 


that is at all close to the desired, the number of antennas must be suffi~ 


cient to prevent ambiguities and give a phase pattern that is geometri-~ 


employed, the steps between antennas must not be greater than approximately 
90°. This requires that the number of antenna be proportional to the 
% aperture when the aperture is very large. When using small apertures, the 


number of antennas must be enough to cause the fundamental of the step phase 














cally consistent. Because of limitations on the phase discriminator 


variation to follow the angle of arrival of incoming signal. The least 
number of antennas that will give a shifting phase variation is three but 


aperture requirements set the lower limit at around 6 or 8 antennas, 
~ 13 - 
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same frequency) is treated mathematically in this report. 


VI ANALYSIS OF EFFECT OF AN INTERFERING WAVE 
(a) Genergl Mathematical Discussion 

When one or more interfering Sinals arrive with random direc- 
tion and time phase with respect to the desired signal, the effect upon 
the Doppler Direction Finder indication is to change the effective 
position of the desired phase variation and thus bring about a bearing 
error. For wide aperture Doppler Systems, the Lurger in amplitude of 
the two signals will make by far the largest contribution to the result~ 
ant phase variation. The stronger signal will "take over" the bearing 
in a manner analogous to the "capture effect! shown by the stronger 
signal in a wide~band frequency modulation system. 

The case which considers a single interfering signal (of the 


A general. 
formula has been developed to give the bearing error (Ao). Fig. (la) 
Plate 6 gives the space diagram illustrating the arrival of the primary 
and interfering waves. The ratio of the interfering signal magnitude 

to that of the primary signal is designated by a, 6 is the arbitrary 
angle of xepsean: of the interfering signal with respect to the other, 

and Y is the time phase difference between the interfering and primary 


Signals, 


For the primary wave 


-@: = — cos( wt + B coswt) (31) 


and for the interfering wave 


€2 = « cos{ut + 8 cos(w,t ~ A) + y] (32) 


cee tee see RSE SRE 





where © = frequency of the carrier 





8 = - electrical aperture of the system 
WW, + rotation or switching rate 
” a ss radius of circle. 
* if e and ey are represented as vectors they can be combined 


as shown in Fig. (16) of Plate 6, The magnitude of the vector, R, is 
of little importance as the Doppler System operates purely on a phase 
coMparison basis. In order to determine the change in phase, the angle 
Ae must be determined: 


a sing 


Ae = are tan ( —————— ) (33) 

L +a cosy 

where o = [wt + B cos(a,t - 9) + y} ~ (wt + B cosw,t] 
= B[ cos(w,t - 9) - cos wt] + (34) 
e The complete expression for the vector resultant of e y and e, 
is: 
R cos[wt + 8 cosw,t + arc tan ( eee (35) 
. 1 + & cosp 

Let yw be the phase angle of this vectors 

yw = wt + B cosw,t + arc tan (ee (36) 
1 +4 cosp 


Without any interfering signal, the phase would be: 
wy = ot + B cos 4t (37) 


Hence, the extra arc tan term must be an error term: 


a sin 9 


Ae = error term = are tan (————— 
; 1 + @ cos@ 


(38) 
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This Ae term will contain some constant terms independent 
of wy (which will cause no phase shift of the B cos w,t term), and 
some cos wt terms as well as harmonics of these terms. The content 
of Ae can be resolved by expanding it in an infinite series as follows 
(see Appendix I for mathematical derivation). 


ad sino 


Ae Tec ) (39) 


are tan ( 


a2 : 
= a sin) - —5~ Sin 29 + 





sin 39 - ... + eae (40) 
The expression for % can be simplified to some extent.. 


Refering to Fig. (1c) 


pg = Bleos (wyt ~ 8) - cosuyt}] + y (34) 
Combining 8 cos(w,t - 9) and -8 cosw,t as vectors gives a vector X 
Fo a 6 iv 
X =v 2=2 cos8 cos(wyt - > -~ a (42) 
a 
pice 8 
=V¥ 2-2 cos 9 sin (at ~—- ) (42) 
Peete eee Q 
Let Bv 2~2cos% = @& and iia) =r 
then 9 = gsinr+ y (43) 
Ae = arc tan [ @ sin(g sin r + ¥) (44) 


1 + a cos(g sin r + y) 





=@ sin(g sinr + y) - sin(Qg sin r + 2y) + see. (45) 


By using the trignometric expansion 


sin(g sin r + y) = sin(g sin r) cosy + cos(g sin r) sin Y (46) 
Ae becomes « sin(g sin r) cosy + @ cos(g sin r) sin y 
2 2 


a 
roe sin(g sin r) cos 2y - or cos(2e sin r) sin 2y +... (47) 


a 


Zp 
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It, is known that 


cos(g sin r) 


Jo(g) + 23 J», (g) cos (2nr) (48) 


sin(g sin r) o zo Jonti( g) sin{(2n + 1)r] 


6 
There is need to consider only the terms involving ( Wot = = ) 


as these are the only terms which can contribute to a phase shift of the 


8 
originel f cos wo, t term. The [2(w,t - aor )] and higher frequency 


terms are generally removed after detection by a band-pass filter, to 
facilitate phase comparison, Let the portion of Ae involving only 
8 
(#1t - -5- ) terms be called Ae! 
i ; 8 
The expansions of cos (g sin r) and sin (g sinr) give (wt - oe ) 


terms only for 1 = 0 


Ae' = af[2J,(g) sin r} cosy + 0 


x 2 


= {2J,(e) sin r} cos2y + 0: 


Prag es eens we hen) pe eO {49 ) 


8 
Met = 2a cosy J, (8 V 2 - 2 cosd ) sin(u,t sre 
2 ee 8 
: -@ cos2y J, (28 v 2 ~ 2 cosO § sin(wit - re, (50) 
8 + 


ee e e e e e 6° e ° ° e oe e e 





ss Iq" 
det = [ 2 ~ (-1)"*? cos ny Js (n8 V2 = 2cos8 Tain (ot ~) (51) 
ona 
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Let the infinite series which is the coefficient of 
sin(w,t - >) be called D. Then 


Ae = D sin(w4t - = ) (52) 


D converges sufficiently rapid. . that only four terms 
( «a, 0, 0°, «*, ) usually need be calculated to give three or four 
place accuracy for D. 
As bearing is given by a shift in the phase of the original 
8B cosw,t term, then bearing error ( Ap ) will result from any uh- 
desired shift in the phase caused by adding the D sin(uyt ~~ ) term. 
The magnitude of this additional shift can be determined by refering to 


Fig. (1d). Let y be the resultant voltage of combining 8 cos mt and 


D Sin (Wt io ore 


y = B cos wt + D sin(wit - (53) ‘A 
= B& cos(w,t — Ap) 6 (54) 
Bearing error = Ap <= are tan [ Denese, (55) 


8-D sin(-A-) 
It should be emphasized that this bearing error assumes a 


linear phase detector and no instrumentation error, In practice a 
phase detector is linear for small phase deviations, but 
approaches a sinusoidal response for larger deviations. 
(b) Figure of Merit for Doppler Systems 
To give some idea of the amount of error encountered in the 
Doppler System under different conditions of the interfering wave, the 
expression (55) was plotted. Three different apertures were chosen 
with the express purpose of showing the effect of increasing the aperture 
and consequent reduction of site errors. In order to compare the 


relative error of the Doppler and Adcock Systems, the expression (23) - of 


Appendix II was plotted under similar circumstantes. (23) is the bearing 


~ 18 - 
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; (5 
error for an Adcock System whose aperture is very small in wavelengths, °° 


The curves are divided into the following groups: 


1. Plates 7, 8, 9,-~Adcock-—Small aperture 
2. Plates 10, 11, 12-~—Doppler—— 4/8 .Aperture 
. 3. Plates 13, 14, 15-—~Doppler—— Aperture 
4. Plates 16, 17, 18—~—Doppler—~5A Aperture 
The word Aperture above, as throughout the report, refers to 
the electrical radius of the system. 
These groups are further divided into different ratios of 
signal magnitudes, namely a = .2,°.5, and .& Different values of 8, 
space separation of the waves are plotted against y , the time phase 
difference. 
A general observation of the curves will reveal the following 
| facts: 
Be 1. Maximum error occurs at y = 0, 180° 
2. At two points the bearing error is zero as was pointed out 
earlier in the report. 
3e As the aperture of the Doppler increases the maximum error 
decreases, 
Comparing the curves for Adcock and Doppler (8 = 45°) it will 
be noted that there is a marked similarity. Since an aperture of 45° 
"is reasonably small, it is just what is expected when the phase front 
curve is considered. It is interesting to let the Doppler Aperture 
approach zero as a limit and then compare the resultant mathematical 
expression with that for the error of the Adcock. From Appendix II we 
extract (31) which is the bearing error of the Doppler as 6B ~ 0 


a sin 8( cosy + a) ( 1+ cosy + a cosO + & cos® cosy ) 
Se ee Pa TT Te ae ar I ee 
4. + Qu cosy (1 + cos9) + 20° cos8 (2 + cos Y cos@) 





Ao= %are tan [{ —— 


] 


+ 20" cosy (2 cos’ 8 + cos - 1) + a* (2 cos*8 - 1) (56) 


1S 
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and from the same appendix the expression for the adcock 


System (23), 


Oe = ¥% are tan [ 


2a sinO(cosy + a@ cos@) 
fe Te a 
1 + 2a cos® cosy + @ (cos 8 - sin 9) 


(57) 


when y= 0° or 360° which is the case for the origin of the 
phase front curve of Plate (5) the two expressions (56) and (57) be~ 


come equal, namely, 


a@ sinS + a sin® cos6 8) 
8 2 = | r eee y (5 / 
a Ap = Yarc tan ( eee ee ee 


If y is placed equal to 180° there results another set of equiv- 
Stent equal formulas. 

However, if ¥ is set°equal to any other angle, the two bearing 
error formulas are not e ual. In the derivation of the adcock formula 
attention was paid to the magnitude of the interference pattern as well 
as the phase. at the 0° 4 180°, and 360° the magnitude is symetrically 
disposed and causes no blur in adcock reading, while at other points the 
blur causes the iainimum reading to be shifted. The difference in actual 
bearing between Adcock and Doppler is very slight and for all practical 
cases it would seem that the result could be extended generally to in- 
clude all small aperture systems, the conclusion being that any. small 
aperture system suffers fromthe same site errors. 

yihen comparing Direction Finding Systems it becomes necessary 
to formulate some standard by which a quick comparison of site error 
reduction can be made. a figure of merit (or demerit) used by H. G. 
Hopkins of the National Physical Laboratory (Hngland) is the following: 
The maximum crror is found for each angle of arrival of the interfering 
signal as the time phase varies from 0° to 360°. The rem.s. value of 


these maximum crrors is used as the figure of demerit of the system, 
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Following this method, (with the aid of Plates 10 through 18), 
the r.mvs.- error angle for the Doppler under different aperture has 


been calculated and appears in the following concise tables. 





TABLE I | 


The important observation of the above table is the large 
reduction of error which is obtainable in an aperture of only one wave- 
length. As was pointed out earlier, the errors for Doppler ( 8 = 45°) 
and Adcock were almost identical, so the above table indicates the 
superiority of wide aperture Doppler over the Adoock in this respecte. 
Theoretically, when the aperture is increased to 8 — 1800° there is a 
further reduction of error, However, there is an upper limit on the 
aperture (due to the constant instrumentation error) above which it be- 
comes unprofitable to further increase 8 .. From the table it would 
appear that éhis limit lies somewhere between 6 ~ 180° (a = d/2) and 
B = 360° (a =). Most of the site reradiation signals have amplitudes 


which are much less than 1 (a0<<1) and this further helps the 
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situation. - - 

It igs interesting to compare Doppler and Adcock systems under 
conditions where the interfering wave is arriving from almost the same 
direction as the desired signal and has comparable magnitude (this is 
the usual case of sky-wave multipath propagation). Apply the following 
‘conditions to expression (57), which is the error for the Adcock with 
stall aperture. Let y= 180° and @-A9 y= 180° is a point of maximum 


error as can be seen from plates (7 through 18). The result is: 





OAD a 
Qe % - Aare tan ( ae ) (59) 
Apply the conditions to expression (55), the error for the 
Doppler. Make the assumption that the aperture of the Doppler is very 


small and the result is: 





Ap x = are tan ( — ) (60) 
Comparison of (59) and (60) shows that for small apertures the srrors of 
the two systems have the same order of magnitude. In particular for 
nearly equal to unity the errors approach 90° in both cases. It is con~ 
cluded that when the apertures of the two systems are small, the bearing 
error for the case of multiple path transmission is the same. However, 


when the aperture of the Doppler is increased the expression for error must 


be approached along different lines. Gonsider expression (55), 


8 
‘ Ao-= are tan pe (55) 
B-D sin) 
2 
" and let 9 = A@ 
Ap x are tan (+) (61) 
a 08% 
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oc Ig" 
where )D = 2 Zz (-1)"*? cos n y¥ dy (n8 V2 - 2 coshd ) 
n= n 
when y = 180° 
© 9qn 





cent Ji (n® ¥2--2 cosh); 


If «<i the summation will converge to a finite value that is 
not directly dependent on ®. ThusAp from (61) can be made to de- 
crease by increasing . The conclusion is that discrimination against 


the weakest of two signals arriving from almost the same angle can be 


made arbitrarly great by increasing the aperture, 
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VIL‘ CONCLUSION 

It is concluded that a Direction Finder operating on the Doppler 
Principle can offer a very practical and accurate solution to the 
problem of determining the direction of arrival of wave energy. shen 
the dimensions of the antenna system are made comparable with that of 
the adcock, it seems to offer no improvement over the later in the 
matter of site error reduction. However, the Doppler System lends itself 
admirably to wide-aperture operation under which conditions it can be 
made to give marked reduction of bearing error caused by site re~- 
radiation. Of course, as a wide aperture system, its use would be 
limited to ground installations and’ shipboard for the medium and high 


frequencies, with the possibility of aircraft application only at very 
and ultra-high frequencics. 
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VIIE .. ‘APFENDIX I 
COMBINATION OF PHASE MODULATE) WAVES 


A more rigorous analysis of the combination of the vactors 


e, = cos (wt + B cos 4,t) 


it 


Go o cos [wt + B cos(w,t - 8) + ¥] 


is presented. 


e, = coswt cos (6 cos Wt) ~ sin wt sin (B cos w,t) 


ep = @ coswt cos{B cos(w,t ~ 9) + vy] - asinwt{sin® eos (w,t-6) + y] 


e, * ep = cosut { cos(% cosw,t)-+a cos{ B cos (wt - 9) + Y 1} 


- sin wt{ sin(B cosa,t) + @ sin[B cos(mt ~ @) + yl} 


In combining two expressions 
BE = A cos wt + B cos wt 


= VATE Be cos (wt + y) 


where 


a 


: 5B 
W = are tan (—~) 


We are mainly interested in the phase of the resultant expression 


( yw). 


in (8 + ; = 
wos ut + axe tan{ tin cosust). + a sin[f cos(ust = 6) + vi, , 


-— cos (8 cosw,t) - @ cos[B cos(w,t - 8) + ¥] 


which can be written 


sin(& cosw,t) + @ sin(6 cosw,t) cos (GB cosw,t) cos[B cos (w,t-6) + y] 
~ o sin {8 cost) cos (8 cosw,t) cos[R cos (wit-8) + y} 
+ sin, (6 cosw,t) sin[8 cost - 9) + y J 
+a cos (8 cosi,t) sin{B cos (wt — 6) + y] 

cos (8 cosw,t) + a@ sin(® cos wt) cos(8 cosw,t) sin[B cos (,t-9) + y] 
~ O sin (6 cosw,t) cos(8 cosw,t) sin[8 cos (w,t-9) + y] 
+a cos, (6 cosw,t) cos[B cos (wit ~ 6) + y]} 
+ @ sin’ (B cosw,t) cos{[B cos (wt ~ 8) + ¥] 


w = Wt + are tan { 


one 


~ a 


Divide the aumerator and denominator by: 


Then 


a 3 


W=wWt+ are tan { 


But 


het 


and 


Then 


cos (8 cosu,t) {1 + & cos(B cosu,t) cos[(B cos(w,t - 6) + y] 


are tan ( 


naa ee || 


+ @ sin(B cosw,t) sin[B cos(w,t — 9) + y}} 


tan(R cosw,t) + 


o tan(B cosw,t) &'sin[B cos(w,t -— 8) + y — B cosu,t] 
1+ cos[B cos(,t — 9) + y - 8 cosw,t] 





x+y ) 
l-xy 


= are tan x + are tan y 





& sin(8 cos(w,t - 8) + ¥ - B cosw,t] 
1 + a cos[ 6 cos(w,t — 8) + y — B cosw,t], 


tan(B cosw,t), 


Bcos (wt - 8) + yY — B cosuyt 


wt + arc tan(tan B cosw,t) + are tan y 
4 


a sin @ 


wt + B cosw,t + arc tan( i 
& cos? 


2 6 - 


a sin{B cos (m,t ~ 9) + ¥ ~ B cosu,t] 
1 +a cos(8 cos (wt — 9) + y — B cosu,t]} 


From identities it is known that 








1 1+ ix 
are tan xX = =~ InG--==s===) 
23 1 - ix 
Applying this: 
1% sing 
" l+ Ltéco 
wv = wt + B cosu,t + ie 
1- 1@ Sing 
1 + a cos@ 
l 1 + @ cosy +2: &sing 
yw = wt + SW,t + ——— |} —— 
oo a ( 1 + a cos - id sing 


3 oie i 
w= ot + 8 cosu,t + — [ In(l + ae™) - In (1 + ae? )j 


21 


No value of @ can make e2? greater than + 1. 


Also q was chosen ~l<q sk . Under these conditions: ae? cannot 


be greater than lor less than -1. 








4 


w= ot + B cost + —— [( cel? Yy7_72? + + oe? +.) 
; 
~(ve? = % oe 70 + Fg? ePIP- + 2. )] 
2 3 
W = wt + B cos w,t + @ sing —- sin2g + 
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It will be noticed that if a = O (no interfering wave), the first 


two turns will give the bearing, However, any interfering signal will 


cause the infinite series to have a values and thus produce error. This 


error is the deviation of phase from the desired as the single element 


progresses around the circle. Analyzing any one of the infinite series 


terms 


é a sin 9 
elimi: 


© aa t cos(w,t - 6) + y - 8 cos o,t] 


ee ey 


& B Wi one Y 


it is found that the error is due to (¢4,8,w, 9, and y ).. The factors 
which carry the greatest weight area,8, and y~. The final bearing 
error will be found by extracting all terms of the expansion containing 
w® yt and it will be found (Appendix II) that the reduction of error is 
mainly dependent on 8 (the aperture of the system). 
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APPENDIX IT 
BEARING ERROR (DUE TO SITE RERADIATION) OF THi DOPPLER SYSTEM AS THE 
APERTURE APPROACHES ZERO (8 - 0), aND COMPARISON WITH ADCOCK ERROR. 


fy s AS previously shown,the general formula for the bearing error of a 


Doppler system is given by: 


| 8 

| 

L Ao = are tan Des J (1) 
8 ~D sin¢®) 

2 

where 

| eo n fy ard 

| Deco se ee ea (2) 
} n=1 n 

I 

i! and 

! J, (ng) = Js BV2-2 cos6 ) (3) 


For any value of n which can be chosen 8 can still be made small 


enough to make valid the following approximation 


Ja (ne) x 2 (4) 


, 
, 
ts 
| 
| 





This causes D to now becomes 


es A r_y\ntl 

n= 23 ios nycyr Sy (5) 
n= n 2 
| 
B= 3 oer ye C1) (6) 
| 
i 
Dy = (aeesy-o ces 2y +o cos BVte. eee. aa (7) 

The series can be synthesised as follows: 
It is true that:* 
; : ; l- a" ee ] = 2 (a” cos2y + a cos4y + a” cos6y + * ev 8 ) 
Se ee 
] as 20° 3 +a 2 

fi sa where a@ < 1 (8) 
‘ 
| 
i * (Dwights Table of Integrals). 
4, & ar 29 ~ 
| 
og Se eee per eG hike a, ae ‘ ey Siig nw cs gee Sad == ws 
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eee ee ee; ~~ a —- eee aes ST Ie ‘ 
(1-4 ) e 2 4 6 
1 ee =~2(a° cos?y-a cos4y-@ cos6by.....) 
where a? < I] (9) 
. 1 ~@ co: 
a = acosy + a? cos2y + a° cos3y +... 6 » 
‘ ‘ : where a? <1 (10) 
adding the last two identities 
- (1-a*) 1 - @ cosy 
re a rege oe ee “ 
eo aes lt+a Sopa: (21) 
= (a cosy - a cos2y + a? GCOBOY teh ee 
Thus: 
a 
a -] 1 ~ a cosy 
D, = e@& (—-—--3--—3 Seis Hp sme ee reat eS 12 
’ g( 1+ 20° +a - 40° cos” Y 1+ a” — 20 cosy ) ( 
The importance of this equation is that it reduces the previously given 
bearing error ( Ap ) expression to an analytic function, 
/ ( 8 ) 
Ap = arc tan [- Ds 2. ] (23) 
8-D, sin( @) 
2 
To simplify by rewrite it: 
4 ? 
a -i ( L- @ cosy) (1+ a+ 20 cosy ) ) 
D, = —— pretenses Sb an ee te agg ee ern en net pan (14 
tae (1 + a”)*— 4a’ cosy (1 +a” — 20 cosy)( 1 +a + 20 cosy ) 
- Combining 
Soe _ +a” + a cosy - 20" cos*y - a” cosy 1 . (15) 
aoe (1 + a)” - 40° cosy 
See aah a ee coe) (16) 
: 5 Cie = cosy)( 1 +a” + 2a cosy ) © 
a + cosy 
, ig ( sea geet (1.7) 
% Bl 1 +a” + 2a cosy 
0, 
Sey oS Desh ae ee) (18) 
] + a@ + 2a cosy 
30. 
ee ea oe ee a oy 
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_ ae a + cosy : (19) 
D, = 208 sin( 5 Nias ae on ca 


(5) 
From the Thesis by R. W. Annis, mentioned in the text, the bearing 
error for the Adcock System is 


oe Pye + Pe (20) 


ne 2n 
If placed in our notation, it would appear 


a sin(@ = y) 


: a sin(@ + y) (21) 
z - Pena OPE) EN Oe eae 1 Me wena tN eae ee he 
One 7 “are tan [ erence yy! + ¥ arc tan[ i amsen 
which is 
: : asin@-y) , %sin(@ + y) 
oe , ite cos(6-~ y) _ 1 + a cos(6 + Y) } (22) 
= rr ee 
ne eae 1 % sin(® - y) sin(@ + y) 
[1 + « cos - y)J {1 + @ cos(6 + y)] 
Expanding and reducing: 
2a sinO@( cosy + & cos6) (23) 
se ed arene eS CS Ee Nee 
ote * Aare tan | 1 + 20 cos® cosy + @ (cos 6 - sin 0) 
Putting the Doppler error expression in this form, the result is: 
8 Boh os 
‘A ¥ ne 8 Dy cos (~~) -~ D,% sind (24) 
= Sea Sa a 
e 2 are wen BR" — 2B D, sin( 6) — Ds cosd 
oe 


Now it is necessary to place the expression for D, (49) in (24) and 
compare (24) with (23). 


Expanding the numerator of (24) 


nn 


l+qa° + 2a cosy 


(25) 
2 ,a% + cosy : 
go 6 
MWephee =. 


8 8 we e0 a + cosy 
D6 cos (-5-) = @872 cos (--) sin (>) ¢ 


Sig hse 


= ated 


ee 


; 9 a” + 2a cosy + cos Y 
Bay ee) ae aie ph ee nel Se 
D, %’sin 4a 8 sin (5) @) stl + 20 cosy) 
4 + cos” 
— 4078" (4% ) (1 — cos6) (sin?) pee ee ee 


(1 + a” + 20 cosy) 


2 2 
eee . a + 20 cosy + cos Y. 
a 8 ( sin® — cos@ sin®) Pie rene (26) 


adding D468 cos ) and —D,” % sin® 
ras «7g? : . 
Numerator = Ried «aaa [(@ + cosy) (1 + a” + 2a cosy) sind 
~ &(sin® — cosO sin) (a7 + 20 cosy + cos’y)] 
0,8” 


5 . 2 
2 be é 
[ fs aes oe a {sin9 (cosy + @) (1 + @ cosy)+ & cos® sinO(a + cosy) "J 


Ly 
2 , 
a 9( cosy + a 
A ee [( 1 + @ cosy + a” cos + a cos® cosy)] (27) 
Expanding the denominator of (24) 
: 9 2 8 a + cosy 
~ 28D —-) = — 28° 20 sin’? (—) ( ——,-—-__—— 
; 8 Dy sin) a sin 3) (45 ae 
a + Gosy 
Se ecosd)> Cae (28) 
a ( cos®) aes eee 
. 2,9 ad + cosy 2 
- p° Qs — 272 sa 2 
cos 4a°8" sin’ ( 5) cos? a oe ats on ee) 
=- 2a8°( 1 - cos®) cos® Ga es (29) 


1 +a” + 2a cosy 
Combining these two terms 


a 8° 
2 = 90) sin(—) = D- 62 [—y———-——.~ 
8 (D sin 5 ) -D" cos 7 Hats Ory) } {1 + 20 cosy(1 + cos6) 


+ 207c0s0(2 + cosy cos9) 
+ 20° cosy (2 cos @ + cos® ~ 1) 
+ a” (2 cos'@ - 1)] (30) 


i BOt 








This gives as the resultant expression for Ap. 


Reh Os tae ee Ser eee: 


Ap dre taal sin®( cosy + &)( 1+ 0 cosy + a cosd + 0 cosé eoey)) (31) 
= d POE ee ae Sie te ee ee eT ee OO oe ee 
° eee ve 1 + 2a cosy( 1+ cos) + 20° cosO( 2 + cos y cos0) * 

+ 20° cosy (2 cos” 4 + cosd — J) 


+ a* (2 cos’ 9- 1) 


Ao was placed in the form (31) in order to compare it with 
equation (23) which comes from the adcock system. It is possible to 
express Ap ina somewhat simpler form by placing the expression for 


D, (19) directly in (1). 


ee! 8 % + cosy 
. eee ——, cart a eae Coane ee mee ee 
208 sin(>-) cos) G75 Dea 
Ap = arc tan{ 5 Fey ] (32) 
RB = 268 sin’ (—) (——»—__—__) 
2 1+ &°+ 2a cosy 


& sin9( a + cosy) 


= tan[ 5 
oes 1+a° + 2a cosy - 2 sin (2) (a + & cos@) ; 
2 


- 


¥ @ sind (a + cosy) (33) 
1 +a cosy + &{% + cosy) cos® 





Ap = arc tanl 
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